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A STUDY  OF  THE  CHARGE  AND  CURRENT  INDUCED  ON 


AN  AIRCRAFT  IN  AN  EMP  SIMULATOR  FACILITY, 


PART  I:  THE  FORMULATION 


A study  of  the  aircraft  skin  current  and  charge  is  made  for  an 
aircraft  under  the  illumination  of  the  electromagnetic  field  produced 
by  a nuclear  EMP  simulator.  This  work  is  motivated  by  the  need  to 
relate  the  response  of  an  aircraft  in  the  EMP  simulator  to  the  response  • 
of  an  aircraft  under  actual  EMP  illumination.  The  analysis  is  directed 
toward  a specific  EMP  simulator  at  Kirtland  AFB,  New  Mexico,  however, 
the  methodology  presented  is  applicable  to  other  types  of  simulator 
facilities. 


r 1 
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I 1.  Introduction 

; An  electromagnetic  pulse  such  as  the  nuclear  EMP  (electromagnetic 

pulse)  incident  on  an  aircraft  will  induce  current  and  charge  densities 
on  its  surface.  These  distributions  represent  the  electromagnetic  fields 
external  to  the  aircraft.  Since  the  interior  fields  are  related  in  a 
linear  fashion  to  the  exterior  fields,  acquiring  a knowledge  of  the 
external  fields  is  the  first  step  in  determining  the  currents  and 
voltages  induced  in  the  electrical  system  of  the  aircraft.  Transfer 
functions,  either  theoretical  or  experimental,  are  used  to  relate  the 
voltage  and/or  current  at  a particular  location  to  the  most  significant 
point  or  points  of  entry  for  the  electromagnetic  energy. 

The  simulation  of  the  nuclear  EMP  requires  the  radiation  of  energy 
over  an  extremely  broad  frequency  band.  At  present  there  are  three 
classes  of  EMP  simulators;  the  bounded  wave  simulator  or  transmission 
line,  the  open  radiator  and  the  hybrid  simulator.^  The  bounded  wave 
simulator  provides  the  best  simulation  of  the  EMP  in  terms  of  frequency 
spectrum  and  field  uniformity.  However,  for  the  illumination  of  very 
large  objects  such  as  aircraft  the  construction  costs  and  the  time 
required  for  construction  become  serious  problems.  Because  of  the  low 
frequency  content  of  the  EMP,  the  dimensions  of  the  open  radiator  type 
of  simulator  become  inordinately  large  and  present  opperational  diffi- 
culties. The  hybrid  simulator  may  be  the  most  cost  effective  of  the 
simulator  facilities.  At  high  frequency  the  hybrid  simulator  behaves 
as  an  open  radiator  and  at  low  frequency  the  near  field  of  the  simulator 
antenna  are  utilized.  However  there  may  be  a tremendous  point-to- 
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point  variation  in  the  field  components  produced  by  this  simulator,  at 
least  to  the  extent  that  the  response  of  an  illuminated  object  may  depend 
upon  the  orientation  of  the  object  as  well  as  its  location  in  the  facility. 

The  Air  Force  Weapons  Laboratory  located  at  the  Kirtland  Air  Force 
Base  uses  various  hybrid  type  simulation  facilities  to  test  the  response 
of  an  aircraft  to  the  nuclear  EMP.^  In  each  of  the  existing  facilities 
the  aircraft  must  rest  on  a concrete  pad  which  eliminates  the  possibility 
of  testing  directly  the  Inflight  mode  response.  Also  the  electromagnetic 
field  components  produced  in  the  facilities  generally  have  a spaclal 
variation  that  is  considerably  different  from  those  associated  with  the 
actual  nuclear  EMP.  Thus  a theoretical  model  is  required  for  the  inter- 
pretation of  test  data  as  well  as  for  the  extrapolation  of  the  test  data 
to  predict  the  inflight  mode  response. 

This  paper  presents  a theoretical-numerical  formulation  for  determin- 
ing the  current  and  charge  distributions  Induced  on  an  aircraft  with  low 
frequency  illumination,  l.e.,  for  wavelengths  that  are  large  in  terms 
of  the  wing  or  fuselage  diameters.  The  analysis  considers  an  electro- 
magnetic model  for  the  aircraft  formed  with  axial  segments  of  bodies  of 
revolution,  where  the  cross  section  radii  and  segment  lengths  are  chosen 
both  for  convenience  and  for  the  most  accurate  physical  representation 
of  the  aircraft  (see  figure  1).  Details  of  the  formulation  are  given  in 
an  earlier  report.^ 

The  foregoing  formulation  yields  only  the  total  axial  current  on 
the  axial  segments  of  the  bodies  of  revolution.  Assuming  a uniform 
circumferential  distribution  can  lead  to  significant  errors.^  With  the 
assumption  of  electrically  thin  segments  the  quasi  state  distribution  of 
current  and  charge  may  be  utilized  to  predict  the  circumferential 


distributions.  This  is  accomplished  using  the  results  from  the 
analysis  of  the  infinitely  long  cylinder. 

In  order  to  determine  the  solution  for  the  aircraft  skin  current 
and  charge  distributions  the  electromagnetic  field  illuminating  the  air-  j 

craft  must  be  known  as  a function  of  frequency.  Experimental  field  map  j 

( 

data  could  be  used.  Two  specific  illuminations  are  considered  incident 
on  the  aircraft.  First  a plane  wave  propagating  normal  to, the  wings  and 
fuselage  is  considered.  Second  the  electromagnetic  field  produced  by 
the  ATHAMAS  (also  called  HPD)  simulator  of  the  Air  Force  Weapons  Labora- 
tory is  considered  (see  figure  2).^’*^  The  later  field  strengths  are 
computed  according  to  the  analysis  of  Blackburn  and  Taylor.^ 

Verification  of  the  theoretical  model  is  accomplished  by  extensive 
comparisons  with  measured  data  that  are  obtained  using  scale  models  of 
aircraft  and  measurements  conducted  on  the  EC-135  aircraft  in  the  ATHAMAS 

simulator.  These  will  be  presented  in  a subsequent  report  (PART  III:  i 

Numerical  Results)  . j 


rv 


2.  Analysis 


a.  Theoretical  Model 


In  solving  for  the  induced  current  and  charge  densities  on  an  air- 
craft the  principal  features  of  the  aircraft — wings,  fuselage,  stabili- 
zers, etc. — are  represented  using  axial  segments  of  bodies  of  revolution. 
For  an  aircraft  such  as  the  EC-135  a convenient  body  of  revolution  model 
is  shown  in  figure  1.  Here  the  bodies  of  revolution  are  either  cylinders 
or  prolate  spheroids.  An  alternative  model  (to  be  discussed  subsequently) 
could  be  constructed  from  segments  of  cylinders.  The  lengths  of  the 
principal  features  of  the  aircraft  are  matched  with  the  model  while 
the  radii  of  the  circular  cross  sections  are  chosen  as  suggested  by 
King.^  The  choice  of  the  radii  is  not  a crucial  factor. 

For  convenience  the  induced  current  and  charge  densities  the  diame- 
ters of  the  principal  elements  or  features  of  the  aircraft  are  considered 
small  in  terms  of  wave  length  so  that  the  axial  current  on  the  elements 
dominates  the  rotational  current  that  may  be  excited.  Therefore,  the 
solution  technique  is  necessarily  limited  to  low  frequency.  To 
determine  the  current  Induced,  the  elements  of  the  aircraft  structure 
are  divided  into  electrically  short  increments  with  an  assumed  sinusoidal 
axial  current,  l.e.,  for  the  nth  segment  of  the  jth  element  with  end 

currents  I.(s.  ) and  I,(s.  ., ), 

j jn'  j j,n+l^’ 

I (g  ) = -J j j jn j ,I}+i i_  (1) 

^ ^ sin  k(s.  ,,  - s.  ) 

j,n+l  jn" 


is  the  induced  axial  current  distribution.  The  unknown  end  currents  of 
the  aforementioned  electrically  short  current  segments  may  be  obtained 
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by  using  the  so-called  extended  boundary  condition  to  develop  a complex 
system  of  linear  equations  and  then  using  standard  numerical  solution 
techniques  to  solve  the  system  of  linear  equat ions . ^ ^ 

The  foregoing  solution  technique  is  numerically  stable  and  is  most 
efficient  for  aircraft-like  dimension.  However,  it  would  be  totally  un- 
suitable for  treating  thin  wire  configurations  due  to  the  convergence 
problems  that  would  occur,® 


b.  Surface  Current  and  Charge  Densities 

In  comparing  the  results  of  the  body  of  revolution  formulation  with 
experimental  data  the  axial  current  that  is  obtained  from  the  formulation 
must  be  related  to  the  Induced  surface  current  density.  Since  the  axial 
£ -igments  of  the  bodies  of  revolution  are  required  to  be  electrically 
thin  it  is  assumed  that  the  circumferential  variation  of  the  current  is 
related  to  the  axial  current  as  on  the  infinitely  long  cylinder. 

Accordingly  consider  the  circumferential  variation  of  the  axial 
surface  current  density  (or  circumferential  magnetic  field).  At  a 
particular  cross  section  of  an  element  of  the  aircraft  structure  the 
magnetic  field  may  be  expanded  in  a Fourier  series  as 


I(s  ) 
2na,(sJ 


(2) 


where 


f(Sj,e) 


n=l 


[A  (s.)cos  n0  + B (s.)sln  n6] 
n'  j'  n'  3' 


(3) 


Evidently  the  function  f(Sj,6)  represents  the  asymmetrical  portion  of 
the  magnetic  field  that  cannot  be  determined  from  the  body  of  revolution 
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formulation  as  discussed  in  the  foregoing. 


In  order  to  obtain  an  approximation  to  the  function  f(Sj,9)  the 
Infinitely  long  cylinderical  scattered  is  considered.  With  an  incident 
plane  wave  propagating  normal  to  the  cylinder  and  with  the  electric 
field  directed  parallel  to  the  cylinder,  the  total  magnetic  field  at  the 
surface  of  the  cylinder  is^ 


hJ°^(Sj,0) 


” j " cos  n6 

^ TT) 

n=l  h''  ''(ka) 
n 


where  I is  the  total  axial  current,  H'  (ka)  is  the  Hankel  function  of 

n 

the  second  kind,  and  is  directed  in  the  6 = 3ti/2  direction.  There- 


fore an  approximation  to  the  asymmetrical  portion  of  the  magnetic  field 
for  a finite  length  cylinder  is 


f (Sj ,9) 


/ 00  j « 

iH  ^ ^ 

n=l  H^^^(ka) 


For  electrically  thin  cylinders,  i.e.,  (ka)^  < < 1,  (5)  yields  the 
magnetostatic  current  distribution^ 


f(Sj,9)  20  • 


where  the  unit  vector  0 is  defined  in  the  right  hand  sense  with  respect 
to  the  direction  of  positive  current.  Using  (6)  in  (2)  yields  the  total 
circumferential  component  of  the  surface  magnetic  field  for  the 
infinitely  long  cylinder  is 


h;°‘(sj,0) 


+ 20  • H 


where  the  corresponding  current  density  is 


*T  f-  n*"  ^ 


J (s.,e)  = hJ°‘^(s.,0)  (8) 

The  axial  component  of  the  magnetic  field  (corresponding  to  a 
circumferential  current)  also  can  be  obtained  from  the  infinitely  long 
cylinder  formulation.  Again  consider  broadside  incidence  but  with  the 
electric  field  now  directed  perpendicular  to  the  axis  of  the  cylinder. 

The  axial  component  of  the  magnetic  field  at  the  surface  of  the  cylinder 
is^ 

h'^°‘^(s.,9)  = -J.(s.,0)  ^ - j2ka  cos  0]  (9) 

z j 6 j z 

provided  (ka)^  <<  1,  for  the  quasistatic  result. 

For  a plane  wave  of  arbitrary  polarization  Incident  from  the  broad- 
side on  an  infinite  cylinder,  the  wave  is  first  decomposed  into  parallel 
and  perpendicular  polarized  components  then  (7)  and  (9)  applied  directly 
for  the  corresponding  polarizations.  These  results  may  also  be  applied 
to  the  aircraft  problem  provided  the  elements  of  the  structure  are 

electrically  thin.  It  should  be  noted  however  that  for  the  aircraft 
“*inc 

problem  H should  also  include  the  magnetic  fields  scattered  from  all 


the  elements  of  the  aircraft.  But  these  generally  can  be  neglected  and 
accurate  results  still  be  obtained.  Therefore,  the  aircraft  surface 
current  distribution  is  approximately,  on  the  jth  element 
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The  modes  of  current  oscillation  associated  with  the  terms  of  (10) 
may  be  readily  identified.  The  term 


2713^  ^ 

represents  an  electric  dipole  oscillation  of  the  axial  current.  The 
term 


20 


j^inc 


represents  a magnetic  dipole  oscillation,  a circulation  of  current  about 
the  element  axis.  For  the  circumferential  current  the  term 

z • H 


also  represents  a circulation  of  current  but  here  it  is  about  the 
circumference  of  the  element  cross  section.  And  finally  the  term 


j z 


hinc  0 


represents  an  electric  dipole  oscillation. 
For  a body  of  revolution  (10)  becomes 


I,(s.)  1 

lihh  * ' “ • " J ' 


[1  - j 2ka  cos  e]6 


(11) 


where  t is  the  unit  vector  tangential  to  the  surface  of  the  body  yet 
perpendicular  to  0.^*^  The  foregoing  is  used  to  predict  the  surface 
current  density  using  the  axial  currents  obtained  from  the  numerical 
solution. 

The  surface  charge  density  is  obtained  via  the  equation  of  con- 


tinuity. It  is 
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P^(Sj.O) 


= i 


til 


dlv  J 

s 


For  a body  of  revolution 


Hence 


wa.(s.)  I ao 
J J ' 


J (s.,0)  + cos  6(s.)  -5 — [a.(s,)t  • J (s.6)] 
s J J dSj  J j s J 


predicts  the  charge  density  associated  with  the  surface  current  density, 
c . Simulator  Field  Computations 

For  an  aircraft  in  a nuclear  EMP  simulation  facility  such  as  the  HPD 
facility  at  the  Klrtland  AFB,  New  Mexico,  a very  complicated  spacial 
variation  of  the  field  incident  on  the  aircraft  occurs.  The  numerical 
solution,  however,  requires  only  that  the  axial  component  of  the  illumi- 
nating electric  field  be  specified  at  a discrete  set  of  points  over  the 
aircraft  configuration.  Hence  the  solution  for  the  Induced  axial  current 
may  be  effected  in  a straight  forward  manner.  Then  the  evaluation  of  the 
surface  current  density  is  accomplished  with  a knowledge  of  the  illumi- 
nating magnetic  field  components  as  discussed  in  the  foregoing. 

One  advantage  of  using  the  body  of  revolution  formulation  is  that 
the  illuminating  field  values  needed  may  be  obtained  experimentally  and/or 
theoretically  since  recently  Blackburn  and  Taylor**  developed  a formulation 
for  calculating  the  electromagnetic  fields  produced  by  a hybrid  type  of 
EMP  simulator  such  as  the  HPD  simulator.  Their  results  compare  quite 
favorably  with  measured  data.  In  any  event  the  theoretical  results  could 
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t 


—m 


be  used  to  supplement  measured  data.  Only  the  time  domain  fields  were 
measured  directly.**  The  Fourier  transform  or  frequency  data  were 
obtained  by  digitizing  the  time  domain  results  and  applying  the  Fourier 
transform  numerically. 

d.  Ground  Plane  Considerations 

For  an  aircraft  in  proximity  to  a ground  plane  the  procedure  for 
determining  the  total  surface  magnetic  field  is  further  complicated  by 
the  interaction  with  the  ground  plane.  Due  to  this  interaction  the  current 
and  charge  densities  tend  to  be  greater  on  the  side  of  the  aircraft  near 
the  ground  plane.  According  to  the  static  solution  for  the  circumferential 
variation  of  the  axial  current  on  an  infinitely  long  cylinder  parallel 
to  a perfectly  conducting  ground,  the  distribution  of  current  normalized 
to  a unit  axial  current  is*^ 

'^1  - la  (s  )/h]2 

I j’  1 + [a^ (Sj)/h]cos  0 


where  h is  the  height  of  the  center  of  the  cross  section  above  the  ground 
plane,  and  0=0  locates  the  section  of  the  cylinder  nearest  the  ground 
plane.  The  incident  magnetic  field  contribution  to  the  current  density 
should  remain  the  same  for  a cylinder  over  the  ground  plane  as  for  the 
cylinder  in  free  space,  except  that  the  ground  reflected  wave  must  be 
included  as  part  of  the  Incident  magnetic  field.  Thus  the  circumferential 
magnetic  field  component  becomes 


I(Sj)f j(Sj,0) 
2 . 


+ 20 


(-inc  ^ jjref  j 


4 • -r 
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where  H is  the  magnetic  field  of  the  ground  reflected  wave.  The 
fcregoing  expression  is  then  used  to  compute  the  axial  current  density 
on  an  aircraft  in  the  ground  alert  mode. 

To  compute  the  corresponding  charge  density  induced  on  an  aircraft 
in  the  ground  alert  mode  the  equation  of  continuity  is  applied  as  mentioned 
in  the  foregoing. 

e.  Time  Domain  Analysis 

The  time  domain  induced  skin  current  and  charge  will  be  calculated 
for  the  aircraft  in  the  HPD  facility  and  for  the  aircraft  inflight 
Illuminated  by  a plane  wave  EMP.  The  body  of  revolution  formulation  that 
is  used  involves  a numerical  solution  for  the  axial  current  over  the 
entire  aircraft  for  a given  excitation  frequency.  To  obtain  the  currents 
under  pulse  excitation,  the  aircraft  currents  are  calculated  for  the 
frequency  band  of  the  incident  pulse  and  then  Fourier  superimposed.  An 
alternative  approach  that  is  more  advantageous  under  certain  conditions 
is  the  use  of  the  singularity  expansion  method®.  Both  procedures  are  to 
be  considered.  The  surface  charge  density  on  the  aircraft  can  be 
obtained  similarly. 

Since  the  excitation  of  the  aircraft  sking  current  is  considerably 
different  for  the  two  cases  considered — simulator  environment  and  the 
inflight  environment — the  current  and  charge  at  certain  locations  will 
be  calculated  and  compared  as  a function  of  frequency.  Then  corresponding 
comparisons  will  be  made  for  the  time  domain  currents  and  charges. 
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1.  Conclusion 

A theoretical-numerical  analysis  has  been  developed  to  determine  the 
low  frequency  skin  current  and  charge  densities  induced  on  an  aircraft 
illuminated  by  an  incident  electromagnetic  field.  The  theoretical  model 
has  been  verified  by  comparison  of  the  results  with  data  that  arc  obtained 
from  scale  model  measurements.  Along  the  principle  features  of  the  air- 
craft the  theoretical  data  is  generally  within  6 db  of  the  measured  data. 

Pulse  excitation  can  also  be  treated  by  use  of  the  Fourier  transform. 
To  Illustrate  the  theoretical  technique  an  aircraft  is  considered  illumi- 
nated by  a unit  step  incident  pulse. 

An  important  application  of  the  analysis,  at  the  present,  is  in  the 
interpretation  of  skin  current  and  charge  data  obtained  in  an  EMP  simulator 


facility. 


JA 


A.  Refcrencfs 


1.  Baum,  C.  E.,  "EMP  Simulators  for  Various  Types  of  Nuclear  EMP 

Environments:  An  Interim  Categorization,"  Sensor  and  Simulation 

Note  151,  Air  Force  Weapons  Laboratory,  Klrtland  AFB,  New  Mexico, 
July  1972. 

2.  Taylor,  C.  C.,  Chen,  K.  T.,  and  Crow,  T.  T.,  "An  Improvement  on 
Wire  Modeling  for  Determining  the  EMP  Interaction  with  Aircraft," 
Air  Force  Weapons  Laboratory  Technical  Report  AFWL-TR-7A-317 , 
Klrtland  AFB,  N.M.,  October  197A. 

3.  Saucer,  M.  I.,  Latham,  R.  W.  and  Varvatis,  A.  D.,  "Relationship 
Between  Total  Currents  and  Surface  Current  Densities  Induced  on 
Aircraft  and  Cylinders,"  Air  Force  Weapons  Laboratory  Interaction 
Note  19A,  Kirtland  AFB,  N.W.,  August  197A. 

A.  Dana,  T.  A.,  Bumgardner,  M.  K. , and  Skogmo,  P.  J.,  "Horizontally 

Polarized  Dipole  Facility  Evaluation,"  Report  No.  AL-IAA,  EG  and  G, 
Inc.,  Albuquerque,  New  Mexico,  January  1975. 

5.  Blackburn,  R.  F.  and  Taylor,  C.  D.,  "On  the  Electromagnetic  Fields 
from  a Hybrid  Type  of  EMP  Simulator,"  Air  Force  Weapons  Laboratory 
Sensor  and  Simulation  Note  211,  Kirtland  AFB,  N.M. , November  1975. 


6.  King,  R.  W.  P.,  Theory  of  Linear  Antennas,  Harvard  University  Press, 
Cambridge,  Mass.,  1957,  Chapter  1,  Section  7,  pp.  lA-20. 

7.  Crow,  T.  T.  Graves,  B.  D.,  and  Taylor,  C.  D. , "The  Singularity 
Expansion  Method  as  Applied  to  Perpendicular  Crossed  Wires,"  IEEE 
Trans.  Ant.  Prop.,  Vol.  AP-23,  pp.  5AO-5A6,  July  1975. 

8.  Butler,  C.  M.  and  Wilton,  D.  R. , "Analysis  of  Various  Numerical 
Techniques  Applied  to  Thln-Wire  Scatterers,"  IEEE  Trans.  Ant.  Prop., 
Vol.  AP-23,  pp.  53A-5AO,  July  1975. 


9.  Harrington,  R.  F.,  Time  Harmonic  Electromagnetic  Fields,  McGraw-Hill, 
New  York,  1961,  p.  233. 

10.  Taylor,  . D.  and  Wilton,  D.  R. , "The  Extended  Boundary  Condition 
Solution  of  the  Dipole  Antenna  of  Revolution,"  IEEE  Trans,  on 
Antennas  and  Propagation,  Vol.  AP-20,  No.  6,  pp.  772-776. 

11.  Taylor,  C.  D.,  "On  the  Circumferential  Current  and  Charge  Distribu- 
tions on  Circular  Cylinders  Near  a Ground  Plane,"  Interaction  Note 
138,  Air  Force  Weapons  Laboratory,  March  1973. 


Figure  Ij  Body  of  Revolution  Model  of  the  EC-135  Aircraft 
Current  Zonea  are  Indicated  with  Dotted  Lines. 


A STUDY  OF  THE  CHARGE  AND  CURRENT  INDUCED  ON 


AN  AIRCRAFT  IN  AN  EMP  SIMULATOR  FACILITY, 
PART  II:  THE  SINGULARITY  EXPANSION  METHOD 


by 

Terry  T.  Crow 
Murall  Kumbaie 
Clayborne  D.  Taylor 


Prepared  for 

Air  Force  Office  of  Scientific  Research 
Bolling  Air  Force  Base 
Washington,  D.  C.  20332 


December  1976 


The  research  reported  In  this  document  was 
sponsored  under  Grant  number  AFOSR-76-2961 


11 


TABLE  OF  CONTENTS 

Page 

Abstract  1 

1.  Introduction  2 j 

2.  Formulation  3 | 

j 

3.  Numerical  Results  4 i 

4.  Conclusion  11  j 

5.  References  12 

i 


! 

i 


A STUDY  OF  THE  CHARGE  AND  CURRENT  INDUCED  ON 
AN  AIRCRAFT  IN  AN  EMP  SIMULATOR  FACILITY, 

PART  II:  THE  SINGULARITY  EXPANSION  METHOD 

For  an  aircraft  under  the  Illumination  of  the  electromagnetic  field 
produced  by  a nuclear  EMP  simulator  the  natural  resonant  frequencies  are 
significantly  affected  by  the  presence  of  the  concrete  pad  supporting  the 
aircraft.  In  order  to  determine  this  effect  the  singularity  expansion 
method  has  been  applied  to  a system  consisting  of  two  perpendicular  crossed 
conducting  thin  cylinders  oriented  parallel  to  a perfectly  conducting 
ground  plane. 


1.  Introduction 


Using  the  singularity  expansion  method  (SEM)  it  is  possible  to 
determine  the  time  domain  scattering  from  a conducting  object  in  terms 
of  a summation  of  damped  sinusoids.^  In  the  evaluation  of  the  SEM 
solution  the  resonant  or  natural  frequencies  of  the  scattering  object 
are  also  obtained.  In  this  report  SEM  is  applied  to  two  perpendicular 
crossed  conducting  thin  cylinders  located  above  a perfectly  conducting 
ground  plane.  This  configuration  may  be  viewed  as  a crude  model  of  an 
aircraft.  And  the  SEM  solution  technique  can  provide  insight  into 
interaction  of  the  aircraft  with  the  ground  plane. 

A few  numerical  results  are  presented.  In  particular  the  natural 
frequencies  are  determined  as  geometrical  parameters  of  the  wires  are 


2.  Formulation 


As  described  in  earlier  work  the  formulation  is  based  on  the  exact 
field  expression  about  a current  filament  with  a sinusoidal  current. 


The  axial  field  component  is** 

^ _ -s|r-ri|/c 


-s|r-r, |/c 


- I'Crp  - 


ir-r2|/c  I 

I ) 


for  the  complex  Laplacion  frequency  s = o + jw,  the  filament  extends  from 
r = rj  to  r = r2,  I'  is  the  axial  derivative  of  the  current  and 


-2  " *'1 


|r2  - rj 


The  radial  component  of  the  field  is 


p . E(r)  = 

4nesp 


r2)I' (r^) 


-s  r-r,  /c 


-s) r-r2 I /c 


- (r-r.)I’(ri)  7=-z=r 

r-r. 


where  is  the  radial  variable  corresponding  to  the  radial  variable  In  a 
cylindrical  coordinate  system  with  the  z axis  coincident  with  the  axis  of 
the  current  filament. 

The  system  of  algebraic  equations  for  the  current  distribution  is 
derived  by  dividing  the  wires  into  segments  and  expressing  the  current  on 
the  kth  segment  of  the  mth  wire  as 


r ^ *<»  ^ » 
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sinh[s(X|^^j^  - x)/c]  + slnh[s(x-Xj^)/c] 

sinh[s(xj^^j^  - 


(3) 


\ 1 ^ 1 


At  this  point  one  must  modify  (1)  and  (2)  which  are  free  space  equations 
to  account  for  contributions  due  to  the  image  currents.  Image  currents  are 
assumed  to  have  the  same  magnitude  as  the  object  currents.  For  object 
currents  parallel  to  the  ground  plane,  image  currents  are  antiparallel  to 
the  object  currents;  for  object  currents  perpendicular  to  the  ground  plane, 
image  currents  are  parallel  to  the  object  currents.  The  scattered  field 
components  about  the  crossed  wires  are  determined  as  functions  of  the  Cij^'s 
by  substituting  (3)  into  the  modified  forms  of  (1)  and  (2).  The  crossed 
wires  are  assumed  to  be  illuminated  by  an  incident  field  whose  electric 
field  is  directed  parallel  to  the  y-axis  (Fig.  1).  Again  one  must  account 
for  the  ground  reflected  wave,  and  the  total  electric  field  illuminating 
the  aircraft  is 

— inc  "Srlnc  :rreflected 

fc.  ^ , — t + t 

total 

“"inc  ' ' Qcsti 

The  total  electric  field  (E  , + E ) is  required  to  satisfy  the 

total 

boundary  condition  that  its  component  tangent  to  the  wire  axis  be  zero. 
Requiring  that  these  equations  be  satisfied  at  a discrete  set  of  match 
points,  one  finds  a linear  set  of  equations  for  the  The  match  points 

are  taken  to  be  the  end  points  of  the  current  segments  referred  to  in 
(3).  Three  comments  are  in  order:  1)  the  King-Wu  charge  boundary  con- 
ditions® replace  three  of  the  above  equations  for  2)  the  Klrchhoff 

current  law  replaces  a fourth  equation;  3)  the  current  at  the  second  match 
point  on  each  segment  is  assumed  to  be  the  average  of  that  on  the  first 
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and  third  match  points — the  first  match  point  on  each  segment  is  at  the 
junction  of  the  four  wires.  The  current  is  assumed  zero  at  the  end  of 
each  wire. 


The  resulting  system  of  linear  equations  for  the  current  can  be 
written  as 

Z(s)  J(s)  = E(s)  (4) 

where  Z(s)  is  the  square  system  matrix,  J(s)  is  a column  matrix  whose 

elements  are  the  a,  's,  and  E(s)  is  a column  matrix  whose  elements  are  the 
k 

incident  field  components  evaluated  at  the  match  points.  Certain  elements 
in  the  E(s)  matrix  are  set  equal  to  zero  as  a result  of  the  boundary  con- 
ditions and  averaging  discussed  in  the  previous  paragraph.  SEM  results 
assert  that  the  current  vector  may  be  written  as 

a a 

E (s)  is  the  Laplace  transform  of  the  incident  waveform,  and  the  s 's  are 
o a 

the  complex  natural  frequencies  defined  as  the  zeros  of  the  determinant 
of  the  system  matrix 

det[Z(s^)]  = 0 (6) 

the  V 's  are  the  natural  modes  defined  by 
a 

Z(s^)  = 0 (7) 

and  the  n^'s  are  the  coupling  coefficients  defined  as^ 

(8) 


J 


the  M 's  are  solutions  of 
a 


Z(s  u = 0 

a a 


where  <a,b>  = 


a • b dt, 


wires 


and 


3Z(s) 

3s 

s=s 


(9) 
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The  details  of  such  calculations  have  been  discussed  in  various  reports.  ’ 
It  is  of  Interest  to  report  that  the  s^'s  tabulated  later  in  this  report 
were  located  using  subroutine  SEARCH^  based  on  the  Cauchy  integral 
theorem  and  verified  using  a Muller  iteration  method.^ 

For  a step  pulse  incident  on  the  wires, 

E 

E (s)  = -^  (10) 

o s 

The  time  domain  current  becomes 

^ s t 

J(t)  = u(t)  Z n — e “ V (11) 

a ® ®(x  “ 


where  u(t)  is  the  unit  Heaviside  step  function. 


3.  Numerical  Results 


r 


Geometrical  parameters  for  the  crossed  cylinder  configuration  are 
illustrated  and  defined  in  figure  1,  A limited  parametrical  study  is 
presented  to  establish  the  feasibility  of  the  method  and  to  determine 
trends.  In  table  1 the  complex  natural  frequencies  are  given  for  crossed 
cylinders  with  aircraft-like  dimensions  in  free  space.  These  form  a 
reference  set  to  compare  with  corresponding  data  for  crossed  cylinders 
loacated  over  a ground  plane.  Data  for  a few  configurations  are  given  in 
Table  2. 

When  S.J  = ilj  the  crossed  cylinder  interaction  is  decoupled,  i.e. 
the  axial  currents  induced  on  the  cylinders  are  the  same  as  would  exist 
on  a single  cylinder.  But  in  general  the  natural  modes  Induced  on  the 
crossed  cylinder  configuration  may  be  separated  into  two  types — coupled 
cylinder  modes  and  Isolated  cylinder  modes.  The  isolated  cylinder 
modes  are  those  that  would  occur  on  an  isolated  cylinder  over  a ground 
plane.  Hence  these  modes  are  independent  of  From  the  data  in 

Table  2,  the  a = 2,6,  •••  modes  are  the  isolated  cylinder  modes  while  all 
the  other  modes  are  the  coupled  cylinder  modes.  But  as  ->■  the 
coupled  modes  vanish. 

Note  that  as  h is  decreased  by  1/2  the  real  parts  of  s^^  decrease 
almost  proportionally  while  the  imaginary  parts  of  s^  change  less  than 
10%.  The  decrease  in  the  real  parts  of  s^  correspond  to  an  increased  Q 
and  the  structure's  decreased  capacity  to  reradlate  electromagnetic 
energy.  If  the  ground  plane  were  lossy  tljen  the  real  parts  of  s^  should 
reflect  the  loss  of  energy  into  the  ground  plane.  Experimental  data 
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Figure  1:  Two  perpendicular  crossed  cylinders  oriented 

parallel  to  a perfect  ground  plane. 
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indicates  that  even  for  a concrete  pad  the  real  parts  of  s^  should  not 
be  as  large  as  the  free  space  values.  This  will  be  discussed  in  Part 
III:  Numerical  Results. 


Table  1: 


Complex  Natural  Frequencies  for  Crossed 
Cylinder  Configurations  in  Free  Space 


+ 

£j  /£  j 


L = 2£.2f  L/a2  ~ 20.,  aj  = a2 , 


0.50) 


a 

1 

2 

3 

4 

5 

6 
7 


-0.2935  + j 2.2922 
-0.4203  + j 2.5082 
-0.3346  + j 3.8392 
-0.6609  + j 6,0656 
-1.0065  + j 8.1563 
-0,9528  + j 8,2939 
-1.1143  + jll.0767 
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Table  2:  Complex  Natural  Frequencies  for  Crossed 


Cylinder  Configurations  Over  a Ground  Plane 


(iil  + 

^ - 

L = 2JH2,  L/aj  = 

20. , a J “ 

a 

= 0, 

25. 

s L/c 
a 

h - O.IL 

i\/li  - 0 

.50, 

h ■ 0. ] 

1 

-0.0130 

+ j 

2.1788 

-0,0202 

+ j 

2.4525 

2 

-0.0426 

+ j 

2.7605 

-0.0426 

+ j 

2.7605 

3 

-0.0864 

+ i 

4.7561 

-0.0470 

+ j 

3.9166 

4 

-0,1872 

+ j 

5.6984 

-0.1775 

+ j 

5.9389 

5 

-0.3299 

+ j 

8.2189 

-0.3015 

+ j 

7.7876 

6 

-0.3289 

+ J 

8.3861 

-0.3289 

+ j 

8.3861 

7 

-0.3167 

+ j 

9.5512 

-0.5359 

+ j 10. 7325 

= 0.25, 

i = 0.2L 

= 0. 

,50, 

h * 0.; 

1 

-0.0342 

+ 

j 

2.1190 

-0.0513 

+ 

J 

2.3606 

2 

-0.0898 

+ 

j 

2.5918 

-0.0898 

+ 

j 

2.5918 

3 

-0.2122 

+ 

j 

4.6285 

-0.1021 

+ 

j 

3.7688 

4 

-0.3591 

+ 

j 

5.3673 

-0.3909 

+ 

j 

5.7399 

5 

-0.7195 

+ 

j 

7.9511 

-0.6691 

+ 

j 

7.5812 

6 

-0.6478 

+ 

j 

8.0338 

-0.6478 

+ 

j 

8.0338 

7 

-0.6609 

+ 

J 

9.3583 

-1,2847 

+ 

jll.0712 
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U.  Conclusion 


With  the  limited  parametric  study  that  is  presented  it  has  been 
shown  that  SEM  applied  to  perpendicular  crossed  cylinders  does  provide 
insight  into  the  aircraft-ground  plane  interaction.  This  interaction  is 
particularly  significant  in  understanding  the  response  of  an  aircraft 
to  the  nuclear  EMP.  However,  the  study  presented  is  not  complete  and  is 
reserved  for  future  investigation. 
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From:  Clayborne  D.  Taylor,  Consultant 

Re:  Preposal  entitled,  "Excitation  of  Modal  Fields  In  Parallel- 

Plate  Transmission  Line" 

1.  The  work  being  proposed  should  be  supported.  Having  been  Involved 

In  strip  line  studies  for  many  years,  the  Electromagnetics  Laboratory 
at  the  University  of  Illinois  is  eminently  qualified  to  do  the  work. 

2.  Obtaining  the  modal  expansion  for  the  electromagnetic  field  within 
a finite  width  parallel-plate  waveguide  is  a very  useful  result. 
However,  the  study  proposed  does  employ  a number  of  judicious 
approximations  of  which  you  should  be  aware.  The  illustration 
below  identifies  the  principle  approximations 


Only  the  TEM  mode  is  assumed  for  this  region 


no  reflections  are  Assume  only  a finite  number  of 

assumed  for  this  discrete  modes  exist  in  the  plate 

interface  region.  The  continuous  spectrum 

is  ignored. 


Also  the  propagation  constants  for  the  discrete  modes  are  obtained 
via  a ray-optical  approximation. 
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The  foregoing  approximations  not  withstanding,  the  proposed  work 
is  a very  meaningful  study.  Moreover  it  can  be  accomplished  with 
the  level  of  effort  that  is  proposed.  It  will  provide  the  most 
accurate  estimate  to  date  for  the  high  order  mode  contribution  to 
the  fields  between  the  parallel  plates. 

It  would  be  quite  helpful  if  the  accuracy  of  the  approximations 
used  in  the  proposed  work  were  estimated  and  included  in  the 
proposal . 
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several  days. 

2.  A meeting  called  by  Col  Wheeler  of  NTO  was  attended  where  the 
executive  summary  of  the  EC-135  assessment  was  discussed. 

3.  Several  days  were  spent  in  the  study  of  the  relationship  between 
individual  wire  currents  and  the  bulk  current  Induced  on  a multi- 
conductor cable. 

4.  A meeting  called  by  ELP  was  attended  where  the  multiconductor 
cable  measurements  being  performed  by  MRC  were  reported.  Their 
work  is  directed  primarily  toward  verifying  the  theoretical  model 
being  developed  by  SAI.  It  was  suggested  that  MRC  also  investigate 
typical  cable  lays  and  routing  for  measurement  purposes. 

5.  Proposals  for  studies  of  simulator-test  object  Interaction  were 
discussed  with  Capt.  Roberts.  Reviews  of  the  proposals  are  in 


preparation. 
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To:  Capt.  L.  Robert's,  ELA 

From:  Clayborne  D.  Taylor,  Consultant 

Re:  Proposal  entitled:  Currents  and  Charges  on  Cylinders  in  a 

Parallel-Plate  Transmission  Line 

1.  Most  of  the  work  being  proposed  should  be  supported.  In  particular, 
the  proposed  measurements  should  be  supported  since  Professor  King 
of  Harvard  University  is  the  world's  authority  on  that  type  of 
research. 

2.  The  measurements  proposed  for  the  crossed  cylinders  in  a parallel- 
plate  transmission  line  facility  should  be  quite  useful  in  pro- 
viding an  additional  means  of  evaluating  test  data  collected  at 
the  EMP  simulator  facilities  utilizing  parallel-plate  transmission 
lines.  In  particular,  the  results  are  necessary  for  interpreting 
test  data  for  very  large  test  objects. 

3.  In  addition  to  the  aforementioned  measurements  analytical  studies 
are  proposed  for  a wire  cross  within  a wire  rhombic  structure  and 
within  a parallel  wire  loop  structure.  It  is  stated  that  these 
studies  may  be  necessary  for  interpreting  the  transmission  line 
measured  data.  But  these  structures  although  amenable  to  analysis 
are  considerably  different  from  the  transmission  line  problem  and 
the  fields  produced  are  very  much  more  complex.  It  is  doubtful 
that  the  proposed  analytic  results  would  be  of  more  than  superficial 


value. 
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4.  It  is  very  interesting  that  the  proposed  research  could  not 

possibly  be  performed  with  the  level  of  effort  that  is  proposed. 

For  example  Profs.  King  & Wu  would  spend  only  3.5  man  months,  ^ 

technicians  would  spend  about  18  man  months,  a secretary  would  ! 

spend  3/4  man  month,  and  a programmer  would  spend  about  1 1/2  J 

weeks.  And  that  is  the  total  proposed  effort.  However,  Prof.  King  j 

would  probably  spend  more  time  than  is  proposed  and  he  probably  | 

j 

J 

will  have  students  available  to  help. 
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A STUDY  OF  THE  CHARGE  AND  CURRENT  INDUCED  ON 


AN  AIRCRAFT  IN  AN  EMP  SIMULATOR  FACILITY, 

PART  III:  NUMERICAL  RESULTS 

ABSTRACT 

A theoretical-numerical  solution  technique  is  used  to  determine 
the  current  and  charge  densities  Induced  on  an  aircraft  by  an  electro- 
magnetic pulse.  To  represent  the  aircraft  the  formulation  uses  an 
electromagnetic  model  formed  with  axial  segments  of  bodies  of  revolution. 
Verification  of  the  theoretical  model  is  accomplished  by  extensive  com- 
parisons with  measured  data  obtained  by  using  scale  models. 

Measured  and  calculated  data  are  also  presented  for  an  EC-135 
aircraft  in  the  ATHAMAS  I (also  called  HPD)  simulator  at  Klrtland  AFB, 

New  Mexico.  In  addition  the  EC-135  is  considered  to  be  illuminated 
Inflight  by  a typical  estimate  of  the  nuclear  EMP.  Comparison  of  the 
inflight  data  with  corresponding  simulator  data  is  made  to  investigate 
the  possibility  of  general  trends. 
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1.  Introduction 


This  report  as  the  sequel  to  parts  I and  II  presents  the  numerical 
results  obtained  using  a theoretical-numerical  formulation  for  determining 
the  current  and  charge  distributions  induced  on  aircraft  an  electromag- 
netic pulse  generated  to  simulate  the  nuclear  EMP  (electromagnetic  pulse) 
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emitted  by  the  detonation  of  a nuclear  weapon.  ’ The  analysis  considers 
an  electromagnetic  model  for  the  aircraft  formed  with  axial  segments  of 
bodies  of  revolution,  where  the  cross  section  radii  and  segment  lengths 
are  chosen  both  for  convenience  and  for  the  most  accurate  physical 
representation  of  the  aircraft.  Details  of  the  formulation  are  given 
in  the  earlier  reports. 

Two  specific  Illuminations  are  considered  incident  on  the  aircraft. 
First  a plane  wave  pulse  propagating  normal  to  the  wings  and  fuselage 
is  considered  for  an  aircraft  inflight.  Second  the  electromagnetic  pulse 
produced  by  the  ATHAMAS  (also  called  HPD)  EMP  simulator  of  the  Air  Force 
Weapons  Laboratory,  Kirtland  AFB  is  considered.  The  later  field  strengths 
are  computed  according  to  the  analysis  of  Blackburn  and  Taylor.^  For 
both  the  inflight  and  simulator  results  two  polarizations  of  the  electric 
field  are  considered  incident  on  the  aircraft. 

An  explicit  contribution  from  the  incident  magnetic  field  must  be 
included  to  obtain  accurate  results  for  the  surface  current  density.  This 
is  accomplished  by  combining  the  axial  current  contribution  to  the 
magnetic  field  with  the  incident  magnetic  field  following  the  result 
for  the  infinitely  long  circular  cylinder.  Both  the  current  and  charge 
density  calculations  are  made  and  presented  for  the  frequency  regime  in 


which  the  analysis  is  valid.  And  Fourier  transform  theory  is  used  to 
obtain  the  corresponding  time  domain  behavior. 
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Verification  of  the  theoretical  model  is  accomplished  by  extensive 
comparisons  with  measured  data.  There  are  comparisons  with  a limited 
amount  of  data  obtained  for  the  EC-135  aircraft  in  the  ATHAMAS  facility. 
But  most  of  the  measured  data  were  obtained  using  scale  models  of  the 
B-1  and  the  EC-135  aircraft  and  a physical  reproduction  of  a cylinder 
model  of  an  aircraft.  The  details  of  the  measurements  are  given  in 
other  reports. In  addition  to  these  comparisons  an  extensive  amount 
of  data  is  presented  for  the  EC-135  inflight  and  in  the  ATHAMAS  facility. 


2.  Numerical  Results 


Before  an  extensive  amount  of  data  Is  presented,  a verification  of 
the  theoretical  model  is  made.  First  comparisons  are  made  with  data 
obtained  from  a physical  reproduction  of  a cylinder  model  of  an  aircraft. 

For  the  calculated  charge  density  data  to  be  presented  only  the 
axial  current  density  is  considered  and  the  equation  of  continuity 

applied.  The  procedure  used  to  obtain  the  surface  current  density  is 

<• 

described  In  Part  I.^ 

a.  Intersecting  Cylinder  Aircraft  Model 

In  figure  1 a model  of  the  B-1  aircraft  is  formed  by  using  inter- 
secting sections  of  wires  or  circular  cylinders.  In  order  to  ascertain 
the  accuracy  of  the  presented  formulation  a scale  model  of  the  cylinder 
model  of  B-1  aircraft  was  constructed  and  the  surface  current  density 
measured  at  various  locations  as  a function  of  frequency  for  topside 
incidence.  If  the  theoretical-numerical  data  and  the  measured  data  were 
precise  for  the  cylinder  (or  wire  model)  then  there  should  be  exact  agree- 
ment. In  figures  2,  3 and  4 measured  and  calculated  data  is  compared. 

Also  shown  on  the  figures  are  data  obtained  from  the  body  of  revolution 
model  (see  figure  5).  Here  the  station  number  on  the  fuselage  is  the 
linear  distance  from  the  nose  of  the  aircraft.  Current  density  measure- 
ments are  made  on  the  shadowed  as  well  as  the  illuminated  side  of  the 
aircraft  model. 

From  figures  2,  3 and  4 the  following  observations  can  be  made. 

(i)  The  current  density  on  the  shadow  side  is  completely  dif- 
ferent from  the  current  density  on  the  illuminated  side  for  frequencies 


above  the  fundamental  resonance  of  the  aircraft. 


(11)  The  shadowing  effect  is  adequately  treated  in  the  analysis. 

(Ill)  The  formulation  yields  accurate  predictions  of  the  fuselage 
current  density  for  frequencies  up  to  about  20  MHz  for  the  B-1  model. 

(Iv)  The  body  of  revolution  model  predicts  lower  peak  currents  and 
higher  resonant  frequencies  than  occur  for  the  wire  model. 

(v)  Because  of  the  substantial  difference  In  the  current  densities 
on  the  wire  and  body  of  revolution  models,  the  body  of  revolution  model 
should  be  a considerable  improvement  on  the  wire  model  of  an  aircraft. 

This  occurs  primarily  because  the  body  of  revolution  model  of  an  aircraft 
is  a better  physical  representation  of  an  aircraft. 

It  has  been  demonstrated  that  the  theoretical-numerical  formulation 
adequately  predicts  the  surface  current  density  induced  on  a crossed 
cylinder  configuration  that  resembles  an  aircraft.  The  next  step  in  the 
verification  is  to  show  that  the  currents  and  charges  induced  on  a body 
of  revolution  model  are  an  approximation  to  those  induced  on  an  actual 
aircraft.  To  that  end  comparisons  are  made  with  data  obtained  from 
measurements  obtained  using  scale  models  of  the  EC-135  and  B-1  aircraft. 

b.  B-1  Scale  Model  Data 

In  figures  6 and  14  the  fuselage  and  wing  current  densities  induced 
on  a B-1  aircraft  are  presented  for  a few  frequencies.  Both  the  wings 
forward  and  wings  swept  configurations  are  shown.  In  the  wings  swept 
configuration  the  B-1  is  almost  a delta  wing  aircraft.  From  these 
figures  the  following  conclusions  are  possible. 

(i)  Although  the  theoretical  model  does  not  include  many  ot  the 
structural  details  of  the  aircraft  the  fuselage  and  wing  current  densities 
can  be  predicted  within  about  6 db  over  virtually  the  entire  aircraft  for 
frequencies  up  to  about  20  MHz. 


(il)  Since  the  distribution  of  current  is  accurate  then  the  charge 
distribution  should  also  be  reasonably  accurate. 

c.  EC-135  Scale  Model  Data 


The  EC-135  aircraft  is  shovm  in  figures  15  and  16,  and  the 
appropriate  body  of  revolution  model  is  shown  in  figure  5.  It  is 
obvious  that  the  physical  appearance  of  this  aircraft  is  more  amenable 
to  modeling  by  the  body  of  revolution  model  than  the  B-1. 

Current  density  and  charge  density  data  are  presented  for  the  EC-135 
aircraft  in  figures  17  through  20.  In  order  to  cover  the  entire  fre- 
quency range  with  the  available  Instrumentation  different  size  scale 
models  were  used.®  These  models  were  various  versions  of  the  Boeing  707 
aircraft  and  did  not  have  exactly  the  same  relative  shapes.  This  accounts, 
in  part,  for  the  differences  in  the  data  acquired  in  overlapping  frequency 
regimes  for  different  scale  models.  On  the  basis  of  the  data  presented  in 
figures  17  through  20  the  following  observations  are  made. 

(1)  The  body  of  revolution  model  accurately  predicts  the  fundamental 
resonances  of  the  aircraft. 

(li)  As  a function  of  frequency  the  predictions  for  current  and  charge 
from  the  body  of  revolution  model  are  within  about  6 db  of  measured  results. 

(Hi)  For  the  fundamental  resonance  the  body  of  revolution  model  pre- 
dicts a lower  resonant  frequency  (less  than  10%  difference)  than  is 
observed. 

(iv)  For  the  higher  order  resonances  the  body  of  revolution  model 
predicts  a higher  resonant  frequency  (about  25%  difference)  than  is 
observed . 

(v)  The  difference  in  current  densities  from  the  illuminated  to 
shadowed  sides  of  the  aircraft  are  accurately  predicted. 
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d.  Ground  Plane  Coupling  Data 

When  an  aircraft  Is  in  the  ground  alert  mode  it  Is  resting  on  a 
concrete  pad  and  it  may  be  grounded  by  means  of  one  or  two  ground  straps 
attached  at  convenient  locations.  In  the  analysis  of  the  aircraft  in 
this  mode  the  ground  straps  are  not  considered  for  convenience.  How- 
ever, the  addition  of  ground  straps  in  the  formulation  is  a minor  modi- 
fication. 

Orienting  a cylinder  parallel  to  a perfectly  conducting  ground  plane 
forms  a two  wire  transmission  line,  the  second  wire  is  the  image  of  the 
cylinder.  This  represents  a very  high  Q configuration.  And  a configu- 
ration of  cylinders  or  cylinder-like  objects  parallel  to  a perfectly 
conducting  ground  plane  must  be  expected  to  exhibit  the  same  behavior. 
Hence  an  aircraft  in  the  ground  alert  mode  over  a perfectly  conducting 
ground  should  exhibit  very  high  Q resonances. 

In  figure  21  the  resonant  frequency  versus  aircraft  height  above  a 
perfectly  conducting  ground  is  shown.  For  the  typical  aircraft  height 
the  variation  is  greatest.  Both  measured  and  computed  data  are  shown 
for  the  intersecting  cylinder  model  of  an  aircraft.  In  figure  22  the 
peak  resonant  current  density  versus  aircraft  height  is  shown.  These 
data  indicate  that  current  density  predictions  are  comparable  in 
accuracy  to  the  inflight  predictions.  However,  the  resonant  frequency 
predictions  will  be  more  in  error  than  corresponding  inflight  predictions. 

At  this  point  the  ground  coupling  effect  on  an  aircraft  for  an 
imperfect  ground  can  only  be  surmised. 

e.  Measured  Data  from  HPD  Simulator  Facility 

For  an  aircraft  in  a nuclear  EMP  simulation  facility  such  as  the 
HPD  facility  at  the  Klrtland  AFB,  New  Mexico,  a very  complicated  spaclal 


variation  of  the  field  Incident  on  the  aircraft  occurs.  The  numerical 


solution,  however,  requires  only  the  axial  component  of  the  illuminating 
electric  field  be  specified  at  a discrete  set  of  points  over  the  air- 
craft configuration.  Hence  the  solution  for  the  induced  axial  current  may 
be  effected  in  a straight  forward  manner.  Then  the  evaluation  of  the 
surface  current  density  is  accomplished  with  a knowledge  of  the  illumi- 
nating magnetic  field  component  perpendicular  to  the  respective  aircraft 
elements . 

The  HPD  facility  at  the  Kirtland  AFB  is  illustrated  in  figure  23. 

Two  aircraft  orientations  are  considered;  first,  the  aircraft  fuselage 
is  located  in  the  plane  of  the  loop  with  the  wing-fuselage  junction 
directly  beneath  the  pulser,  and  second  the  aircraft  fuselage  is  oriented 
perpendicular  to  the  plane  of  the  loop  with  the  wing-fuselage  junction 
directly  beneath  the  pulser. 

In  figures  2A  through  27  the  time  domain  axial  current  densities 
measured  at  a few  locations  on  the  aircraft  are  shown  for  the  EC-135  in 
the  aforementioned  parallel  orientation.  Also  shown  are  calculated 
results  obtained  by  considering  two  separate  approximations.  First  the 
ground  beneath  the  aircraft  is  assumed  to  be  perfectly  conducting  (with 
image)  and  second  the  ground  beneath  the  aircraft  is  assumed  to  be  free 
space  (without  image).  The  illuminating  field— with  ground  bounce — is 
considered  to  be  the  same  in  both  cases.  At  the  HPD  facility  the  air- 
craft is  supported  by  a concrete  pad  about  3 feet  thick  with  steel 
reinforcement.  The  constitutive  electrical  properties  of  the  concrete 
are  estimated  to  be,®  conductivity  ~ .008U/m  and  dielectric  constant 
-20  at  5 MHz.  Thus  it  is  expected  that  the  measured  data  fall  somewhere 
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between  (within  ~ 6 db)  the  calculated  data  for  the  perfect  ground  and 
the  no-ground  approximations.  A more  satisfactory  theoretical  treatment 
including  a lossy  ground  is  currently  being  pursued  by  the  authors. 


From  figures  24  through  27  three  principle  observations  are  readily 
made.  First,  the  measured  peaks  are  roughly  a factor  of  two  greater  than 
the  calculated  peaks.  Second,  for  the  very  early  time  behavior  the 
perfect  ground  and  no-ground  data  essentially  coincide  for  the  points  on 
top  of  the  fuselage.  And  third,  the  late  time  oscillation  of  the  measured 
currents  are  more  attenuated  than  the  perfect  ground  results  and  less 
attenuated  than  the  no-ground  results. 

The  Fourier  transforms  of  the  time  domain  measurements  are  presented 
in  figures  28  through  31.  Theoretical  data  for  the  perfect  ground  and 
no-ground  situations  are  also  Included.  In  general  the  magnitude  of  the 
Fourier  transforms  of  the  measured  data  are  less  than  the  calculated 
transforms  by  about  a factor  of  two.  This  is  due  in  part  to  the  inac- 
curacy of  the  Blackburn-Taylor  predition  of  the  simulator  field  components. 
Also  it  is  noted  that  the  no-grcund  results  are  more  similar  in  shape  to 
the  Fourier  transform  of  the  measured  data.  But  most  important  the 
calculated  resonant  peaks  occur  at  frequencies  lower  than  corresponding 
resonances  in  the  measured  data.  This  also  has  been  observed  in  the 
comparisons  with  scale  model  data. 

f.  Computed  Results  for  HPD  Simulator  Facility 

Comparisons  with  measured  data  in  the  preceding  section  Indicate 
the  accuracy  of  the  calculated  results.  In  this  section  the  axial 
current  densities  at  a number  of  locations  over  the  aircraft  structure 
is  presented.  Figures  22  through  35  exhibit  the  simulator  magnetic 
field  at  a few  points  on  the  aircraft  structure  in  a perpendicular 
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orientation  (see  the  figures  for  an  Illustration).  In  figures  36  through 
43  the  magnitude  of  the  Fourier  transform  of  the  axial  current  density 
Is  presented.  The  point  at  which  the  current  Is  evaluated  Is  Illustrated 
by  the  aircraft  figure.  Both  the  perfectly  conducting  ground  (with  Image) 
and  the  no-ground  (without  Image)  results  are  presented.  Figures  44 
through  47  exhibit  the  time  domain  results  corresponding  to  the  afore- 
mentioned Fourier  transform  data. 

Charge  density  data  for  the  EC-135  In  the  HPD  simulator  Is  pre- 
sented In  figures  48  through  56.  Here,  as  before,  the  aircraft  Is  In  the 
perpendicular  orientation.  It  should  be  mentioned  that  for  the  aircraft 
In  the  perpendicular  orientation  there  Is  no  axial  current  Induced  on 
the  fuselage  and  the  vertical  stabilizer.  That  Is  the  reason  for  no 
current  and  charge  density  data  being  presented  for  these  elements. 

In  preparation  for  the  presentation  of  current  and  charge  density  data 
on  the  aircraft  in  the  parallel  orientation, the  HPD  magnetic  field  com- 
ponents are  shown  for  a few  locations  along  the  aircraft  structure. 

Figures  61  through  93  present  current  density  data  at  points  over  the 
aircraft  structure  for  both  the  perfect  ground  and  no-ground  approxima- 
tions. Again  the  position  of  the  current  evaluation  Is  Indicated  in  the 
aircraft  figure.  Figures  94  through  129  present  corresponding  charge 
density  data. 

g.  Computed  Results  for  Inflight  Mode 

As  mentioned  In  the  foregoing  the  current  and  charge  density  induced 
on  an  aircraft  in  a simulator  facility  must  be  related  to  the  corre- 
sponding densities  on  an  aircraft  Inflight.  To  that  end  the  theoretical 
model  Is  considered  to  be  Illuminated  by  the  typical  estimate  of  the 


nuclear  EMP  when  the  aircraft  is  in  the  inflight  mode.  Accordingly  the 
nuclear  EMP  is  approximated  by  a plane  wave  pulse  with  electric  field,’ 
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a = 1.5  X 10®  s"‘ 
e = 2.6  y 10®  s”‘ 


And  the  accompaning  magnetic  field  is 

H(t)  = ^ E(t) 

where  ^ = 120tt  ohms. 

First  current  and  charge  density  data  are  presented  for  the  electric 
field  polarized  parallel  to  the  fuselage  (corresponding  to  the  parallel 
orientation  of  the  aircraft  in  the  HPD  simulator  facility).  These  data 
are  presented  in  figures  130  through  175.  Second,  the  current  and 
charge  density  are  presented  for  the  electric  field  polarized  perpendicu- 
lar orientation  of  the  aircraft  in  the  HPD  simulator  facility.  As  before 
the  location  of  the  current  and  charge  evaluation  is  illustrated  in  the 
figure.  Also  the  Illumination  Is  from  the  topside  of  the  aircraft. 

h.  Extrapolation  to  Free  Space 

The  peak  electric  and  magnetic  fields  for  the  Inflight  mode  illumi- 
nation is  roughly  a factor  of  3 greater  than  the  peak  field  strengths 
produced  by  HPD.  Hence  the  peak  currents  and  charge  densities  are 
expected  to  differ  by  corresponding  amounts.  For  example,  on  the  for- 
ward fuselage  section  at  the  wing-fuselage  junction  the  peak  current 
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density  on  top  of  the  fuselage  in  the  simulator  is  •“  150  A/m  at 
t ~ 0.1  ps  from  figure  89.  And  from  figure  146  the  corresponding  peak 
current  density  for  the  inflight  mode  is  ~ 640  A/m  at  t ~ 0.05  ps.  For 
another  example,  the  current  density  midway  and  on  top  of  the  wing  is 
from  figure  91  is  ^ 60  A/m  at  t ~ 0.1  ps.  And  from  figure  151  the 
corresponding  peak  current  density  is  ~ 280  A/m  at  t ~ 0.025  ps. 

A comparison  of  the  charge  densities  for  the  aircraft  Inflight  and 
in  the  HPD  simulator  is  made  in  a similar  fashion  to  that  for  the  current 
densities.  From  figure  125  the  peak  charge  density  at  the  wing  tip  for 
an  aircraft  in  the  HPD  facility  is  ~ .04  x 10  c/m^  at  t ~ .4  ps.  And 
the  corresponding  charge  density  for  the  inflight  response  from  figure 
172  is  ^ .13  X 10  c/m^  at  t ~ .4  ps. 
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3.  Conclusion 

A theoretical-numerical  analysis  has  been  developed  to  determine 
the  low  frequency  skin  current  and  charge  densities  Induced  on  an  air- 
craft illuminated  by  an  incident  electromagnetic  field.  The  theoretical 
model  has  been  verified  by  comparison  of  the  results  with  an  extensive 
amount  of  measured  data.  In  general,  along  the  principle  features  of 
the  aircraft  the  theoretical  data  is  within  6 db  of  the  measured 
data . 

An  important  application  of  the  analysis,  at  the  present,  is  in  the 
interpretation  of  skin  current  and  charge  data  obtained  in  an  EMP 
simulator  facility.  Typical  results  are  presented  for  the  HPD  facility 
at  the  Kirtland  AFB,  New  Mexico. 


REFERENCES 


1.  Taylor,  C.  D.,  Chen,  K.  T.  and  Crow,  T.  T.,  "A  Study  of  the  Charge 
and  Current  Induced  on  an  Aircraft  in  an  EMP  Simulator  Facility, 

Part  I:  The  Formulation,"  Final  Report  Grant  No.  AFOSR-76-2961 , 

submitted  to  Air  Force  Office  of  Scientific  Research,  December 
1976. 

2.  Crow,  T.  T.,  Kumbale,  M.  and  Taylor,  C.  D.,  "A  Study  of  the 

Charge  and  Current  Induced  on  an  Aircraft  in  an  EMP  Simulator 
Facility,  Part  II:  The  Singularity  Expansion  Method,"  Final 

Report  Grant  No.  AFOSR-76-2961,  submitted  to  Air  Force  Office  of 
Scientific  Research,  December  1976. 

3.  Blackburn,  R.  F.  and  Taylor,  C.  D. , "On  the  Electromagnetic  Fields 
from  a Hybrid  Type  of  EMP  Simulator,  "Air  Force  Weapons  Laboratory 
Sensor  and  Simulation  Note  211,  Klrtland  AFB,  N.  M. , November  1975. 

A.  Knott,  E.  F.,  "Surface  Field  Measurements,  "Air  Force  Weapons 
Laboratory  Interaction  Application  Memo  5,  Kirtland  AFB,  N.M. , 
September  1974. 

5.  Llepa,  V.  V.,  "Sweep  Frequency  Surface  Field  Measurements,"  Air 
Force  Weapons  Laboratory  Technical  Report,  AFWL-TR-75-217 , 

Kirtland  AFB,  N.M. , July  1975. 

6.  Liepa,  V.  V.,  and  Dusette,  D.L.,  "Free  Space  EC-135  Scale  Model 
Measurements,"  to  be  published  as  an  Air  Force  Weapons  Laboratory 
Technical  Report. 

7.  Dana,  T.  A.,  Bumgardner,  M.  K. , and  Skogmo,  J.  P.,  "Horizontally 
Polarized  Dipole  Facility  Evaluation,  "Report  No.  AL-144,  EG  and  G, 
Inc.,  Albuquerque,  N.M. , January  1975. 

8.  Castillo,  J.  P. , "Effects  of  Wave  Reflection  on  Objects  Near  a Plane 
Ground,"  Air  Force  Weapons  Laboratory  ATHAMAS  MEMO  No.  8,  Kirtland 
AFB,  N.  M.,  May  1975. 

9.  Barnes,  P.  R.,  "The  Effects  of  Electromagnetic  Pulse  (EMP)  on  State 
and  Local  Radio  Communications,"  Oak  Ridge  National  Laboratory 
Report  No.  ORNL-4873,  Oak  Ridge,  TN,  October  1973. 


14 


» T • 'V  *• 


Figure  2-  Current  Density  on  the  Topside  of  the  Fuselage  at 
Station  No.  5.0  meters  for  Topside  Illumination 
with  Parallel  to  the  Fuselage. 
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FREQUENCY  (MHz) 

Figure  3:  Current  Density  on  the  Topside  of  the  Fuselage  at  Atation 

No.  22.5  meters  for  Topside  Illumination  with  Parallel 

to  the  Fuselage. 


Figure  4:  Current  Density  on  the  Bottom  of  the  Fuselage  at 

Station  No.  22.5  meters  for  Topside  Illumination 

with  Parallel  to  the  Fuselage. 
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Figure  7:  Comparison  of  B-1  Scale  Mode  Measurements  with 

Theoretical  Results  for  Top  Side  Illumination 
and  Wings  Forward,  f =6.4  MHz. 


Theoretical  Results  for  Top  Side  Illisalnatlon 
and  Wings  Forward,  f » 10.4  MHz. 


Figure  10:  Comparison  of  B-1  Scale  Model  Measurements  with 

Theoretical  Results  for  Top  Side  Illumination 
and  Wings  Forward.  f = 18.16  Mllz. 
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Figure  lA:  Comparison  of  B-1  Scale  Model  Measurements  with 

Theoretical  Results  for  Top  Side  Illumination 
and  Wings  Swept.  f = 18.16  MHz. 


Figure  16:  EC-135  Aircraft  in  the  Inflight  Mode.  Wing  and 
Horizontal  Stabilizer  Station  Numbers  are  Shown 

in  Meters . 
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the  Top  Side  of  the  EC-135  Fuselage 
meters  for  Top  Side  Illumination 
to  the  Fuselage. 


Figure  20:  Charge  Density  on  the  Nose  of  the  EC-135  Aircraft 

for  Top  Side  Illumination  with  Parallel 

to  the  Fuselage. 


(MH 


d (in  m)  / Sp 


Figure  21:  Resonant  frequency  vs.  ground  spacing  for  the 

circular  cylinder  nwdel  of  an  aircraft 
(dimensions  appropriate  to  the  B-1). 
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Figure  22:  Surface  current  density  at  resonance  vs. 

ground  spacing  for  the  circular  cylinder 
model  of  an  aircraft  with  dimensions 
appropriate  to  the  B-1. 
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TIME  (microseconds) 

Figure  25:  Time  history  of  the  axial  current  density  on  the 

bottom  of  the  fuselage  (13.5  meters  from  the 
nose)  for  the  EC-135  aircraft  in  the  HPD  simulator 
facility. 
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Figure  27:  Time  history  of  the  axial  current  density  on  the 

top  of  the  fuselage  (10  meters  from  the  nose)  for 
the  EC-135  aircraft  in  the  HPD  simulator  facility. 
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RIRCRRFT  F8R  TBPSIDE  INCIDENCE  DBUBLE  EXPBNENTIRL  PLRNE  HAVE 


CUAAENT  0EN3ITT  (INCLUDING  HINC  FIELD)  AT  PflSITION  27  ON  THE  EC-I3S 
RinCRHFT  FRR  T8PSIDE  INCIDENCE  DflUBLE  EXPONENTIAL  PLANE  HAVE 
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Figure  150: 


CURRENT  OENSITT  (INCLUDING  HINC  FIELD)  AT  PflSITISN  7 ON  THE  EC-135 
AIRCRAFT  FDR  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  PLANE  HAVE 


CHARGE  DENSITY  RT  PSSITiaN  69  AN  THE  EC-13S  AIRCRAFT 
FSR  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  PLANE  HAVE 


Figure  162: 


CHARCc  DENSITY  RT  POSITION  71  ON  THE  EC-135  AIRCRAFT 


FOR  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  PLANE  HAVE 


Figure  163: 
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CHARGE  OENSITT  AT  PGSITISN  71  ON  THE  EC-I3B  AIRCRAFT 
FOR  TflPSlOE  INCIOENCE  00U8LE  EXPONENTIAL  PLANE  NAVE 
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Figure  165: 


Figure  166: 


CHARCt  DEN3ITT  AT  P6SITISN  SU  SN  THE  EC-13S  AIRCRAFT 
F9R  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  PLANE  HAVE 
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Figure  174: 

CHARGE  OENSITT  AT  PeSITION  7 SN  THE  EC-135  AIRCRAFT 
F6R  TOPSIDE  INCIDEHCE  DOUBLE  EFPONEHTIAL  PLANE  HAVE 
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Figure  175: 
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CURRENT  DENSITY  (INCLUDING  HINC  FIELD)  RT  P8SITI8N  27  8N  THE  EC-135 
RIRCRRFT  F8R  T8PSI0E  INCIDENCE  D8U8LE  EXP8NENTIRL  PLANE  HAVE 
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Figure  178: 

CURRENT  DENSITY  (INCLUDING  HINC  FIELD)  AT  P8SITI8N  30  8N  THE  EC-135 
AIRCRAFT  F8R  TflPSIDE  INCIDENCE  DBUBLE  EKP8NENTIAL  PLANE  NAYE 


T8P  SIDE  8F  AIRCRAFT  HITH8UT  IHAGE 
B8TT8H  SIDE  8F  AIRCRAFT  HITH8UT  IMAGE 
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Figure  179: 
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CURRENT  DENSITY  (INCLUDING  HINC  FIELD)  AT  P83ITI9N  7 SN  THE  EC-13S 
RIRCRRFT  F9R  T9PSIDE  INCIDENCE  D9UBLE  EXP9NENTIRL  PLANE  HAVE 


rHRfiOE  DENSITY  BT  P8SITI9H  13  ON  IHE  EC-135  fllRCRBFI 
F9R  T9RSI0E  INCIDENCE  D9UBLE  E»P9NENnflL  PLRNE  NAVE 
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Figure  184: 


CHRRCE  DENSITY  BT  P9SITI9N  7 8N  THE  EC-135  fllRCRBFT 
FOR  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  PLANE  HAVE 
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Figure  185 


CHORGE  DENSITY  RT  RBSITIDN  33  AN  THE  EC-13' 
FOR  T8PSI0E  INCIDENCE  OflUBLE  EKBflNENtlRL  P' 
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Figure  186: 


CHARGE  DENSITY  AT  POSITION  33  ON  THE  EC-13 
FOR  TOPSIDE  INCIDENCE  DOUBLE  EXPONENTIAL  P 
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Figure  189: 
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Part  II  - For  an  aircraft  under  the  illumination  of  the  electromagnetic  field 
produced  by  a nuclear  EMP  simulator  the  natural  resonant  frequencies  are 
significantly  affected  by  the  presence  of  the  concrete  pad  supporting  the 
aircraft.  In  order  to  determine  this  effect  the  singularity  expansion 
method  has  been  applied  to  a system  consisting  of  two  perpendicular  crossed 
conducting  thin  cylinders  oriented  parallel  to  a perfectly  conducting 
ground  plane. 

Part  III  - A theoretical-numerical  solution  technique  is  used  to  determine 

the  current  and  charge  densities  induced  on  an  aircraft  by  an  electromagnetic 
pulse.  To  represent  the  aircraft  the  formulation  uses  an  electromagnetic 
model  formed  with  axial  segments  of  bodies  of  revolution.  Verification  of  the 
theoretical  model  is  accomplished  by  extensive  comparisons  with  measured  data 
obtained  by  using  scale  models.  Measured  and  calculated  data  are  also 
presented  for  an  EC-135  aircraft  in  the  ATHAMAS  I (also  called  HPD) 
simulator  at  Kirtland  AFB,  New  Mexico.  In  addition  the  EC-135  is  considered 
to  be  illuminated  inflight  by  a typdical  estimate  of  the  nuclear  EMP. 
Comparison  of  the  inflight  data  with  corresponding  simulator  data  is  made  to 
investigate  the  possibility  of  general  trends. 
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